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Abstract 

We treat the noise-activated escape from a one-dimensional potential well of 
an overdamped particle, to which a periodic force of fixed frequency is applied. 
We determine the boundary layer behavior, and the physically relevant length 
scales, near the oscillating well top. We show how stochastic behavior near 
the well top generalizes the behavior first determined by Kramers, in the case 
without forcing. Both the case when the forcing dies away in the weak noise 
limit, and the case when it does not, are examined. We also discuss the 
relevance of various scaling regimes to recent optical trap experiments. 

PACS numbers: 05.40.-a, 05.70.Fh 
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The phenomenon of weak white noise inducing escape from a one-dimensional poten- 
tial well was studied by Kramers |]. If e denotes the noise strength (e.g., e oc k B T in 
thermal systems), and AE measures the depth of the well, then the escape rate A falls off 
like exp(— AE/e) as e — > 0. The case when the trapped particle is overdamped is easiest 
to analyse. If the particle, after each escape, is reinjected at the bottom of the well, and 
a steady state has been set up, then in the interior of the well its position will have a 
Maxwell-Boltzmann distribution. Kramers determined that this distribution must be mod- 
ified near the well top, by being multiplied by a 'boundary layer function' that incorporates 
outgoing boundary conditions. From this modified distribution, he was able to determine 
the weak-noise limit of the escape rate, including the all-important pre-exponential factor, 
by computing the probability flux over the well top. 

The Kramers formula and its multidimensional generalization have been extended in 
many ways 0|§. There have been extensions to non-overdamped particles and to colored 
noise. There have also been extensions to the case when even though the particle is over- 
damped and the noise is white, the noise-perturbed dynamics of the particle fail to satisfy 
detailed balance. This may be due to localized 'hot spots' |4j] or, in multidimensional sys- 
tems, to nonconservative deterministic dynamics 0. 

However, there is one experimentally important case that has not been exhaustively 
studied. That is when the system parameters are periodically modulated. A full analysis 
of escape driven by weak noise, in such systems, would shed light on the Kramers limit 
of stochastic resonance. It would also clarify the effects of barrier modulation on phase- 
transition phenomena. 

It is now possible to construct a physical system (a mesoscopic dielectric particle that 
moves, in an overdamped way, within a dual optical trap 0) that provides a clean exper- 
imental test of the three-dimensional Kramers formula. The rate at which thermal noise 
induces escape agrees well with the predictions of the formula. Adding an external force, of 
fixed period rp, would yield a periodically modulated system M, of the sort that has not 
yet been fully analysed. A complete treatment of escape from a well of a periodically driven 
overdamped particle, or equivalently the escape of an overdamped particle from a 'sloshing 
potential well', would be desirable. 

Smelyanskiy, Dykman, and Golding treated this phenomenon perturbatively, in one di- 
mension ||. They derived a Kramers prefactor incorporating /, the strength of the periodic 
forcing. It applies if the ratio f/e is set to a constant as e — > 0. That is, the forcing is taken 
to die away in the weak-noise limit. Lehmann, Reimann, and Hanggi || treated nonpertur- 
batively the case when / is independent of e, using path integral techniques, and worked out 
a numerical scheme for computing the /-dependent prefactor. They also examined the 'in- 
stantaneous escape rate', which in the steady state is a r^-periodic function of time. In a 
simulation of a special case (a well with a perfectly harmonic top), they noted that in the 
weak-noise limit, the maximum of the instantaneous escape rate cycles slowly around the 
interval [0, t f ). 

In this Letter, we go beyond M and ||. By treating the case / oc e a , where a is an 
arbitrary nonnegative power of e, we determine the relation between their respective scaling 
regimes. In the weak- noise, weak-forcing limit, there are three physically important length 
scales near the oscillating well top, of sizes proportional to e 1//2 , /, and / 1//2 . Crossover 
behavior will result if / oc e 1 / 2 , and the case / oc e can itself be viewed as a crossover regime. 
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When / is independent of e, we use facts on noise-induced transport through unstable 



limit cycles to illuminate the 'cycling' phenomenon [10]. At any t in [0,t>), the normalized 
instantaneous escape rate oscillates periodically in loge as e — > 0. We supply a formula for 
the period, and give a physical explanation for the logarithmic slowness. 

More importantly, we place the case of e-independent periodic forcing firmly in the 
Kramers framework, by determining how the Maxwell-Boltzmann distribution is modified, 
in the boundary layer of width 0{e 1 / 2 ) near the oscillating well top. As / — > 0, it approaches 
the modified distribution of Kramers 0. The case when / oc e in the weak-noise limit is 
intermediate between the case of e-independent forcing and the case of zero forcing, and its 
boundary layer behavior is intermediate too. 

Scaling Regimes. — Initially, we work in terms of dimensional quantities. The Langevin 
equation for a driven Brownian particle in a potential well U = U(x) is 



mi + ■ymx = —U'(x) + Fv(t) + \l2m r ykBT rj(t). (1) 

Here 7 is the damping, F a dimensional measure of the driving, v a dimensionless periodic 
function of unit amplitude, and 77 a standard white noise. In the overdamped (large-7) limit, 
the inertial term can be dropped, leaving 

x = - V \x)+fv(t) + ^r l (t). (2) 

Here V = U/jm, f = F/jm, and e = 2ksT / r ym. 
The Kramers formula for the / = escape rate is 

X ~^exp(-AU/k B T) 

JV"(x 8 ) \V"(x u )\ 
= V , exp(— A£//e), (3) 

Z7T 

where uj s = JU"(x s )/m and u u = yJ\U"(x u )\ jra are the oscillation frequencies about 
the bottom x s and top x u of the well, and AE = 2AV. Eq. (|3|) follows from 
Kramers's modification of the steady-state Maxwell-Boltzmann weighting exp [—U(x)/kBT], 
i.e., exp [— 2V(x)/e]. If n denotes the inward offset from x u , his modifying factor is 

erfc[-n/y / e7jy^K)l]. 

If / 7^ 0, there are two regimes, depending on the size of / as e — > 0. Since [e] = [AE] = 
L 2 /t and [/] = L/t, where L denotes length and t denotes time, comparing / with e must 
be done with care. / will be 'small' or 'large' in the Kramers limit if it is small or large 



compared to a quantity with dimensions L/t, namely J\V"(x u )\ e. 

In physical terms, there are two regimes because there are two length scales at the well 
top, and one or the other is larger. The first is the length scale in Kramers's modification. 
There is a boundary layer of width w ^2ksT / \U"(x u )\, i.e., yJe/\V"(x u )\, within which 
'physics occurs'. This (9(e 1//2 ) quantity is the diffusion length: the distance from the top to 
which the particle must approach, to acquire a substantial chance of leaving the well rather 
than falling back. 

If a periodic force is applied, a second length scale becomes important. The top of 
the well will oscillate periodically around the unperturbed top x u by an amount roughly 
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equal to Fj \U"(x u )\. If this length scale is substantially smaller than the first, to a first 
approximation the boundary layer will not oscillate. But if the opposite is true, escape 
dynamics should be strongly affected by boundary layer oscillations. The crossover occurs 
when F ~ \Jl \U"(x u ) \ ksT, i.e., when F ~ \jlmijj\kBT . In normalized units, this criterion 

is / « ^\V"{x u )\e. 

So if / oc e a in the Kramers limit (e — > 0), a > 1/2 and < a < 1/2 belong to different 
regimes. The a = 1 results of Ref. || presumably extend to the entire a > 1/2 regime. 
Similarly, our treatment of a = below could be extended to cover the < a < 1/2 regime. 
These two regimes should be kept in mind when conducting experiments on noise-driven 
escape in periodically driven systems. In the Kramers limit, only when the forcing F is 
much less than y^mu^/c^T is a simple perturbative modification of the Kramers formula 
likely to apply. 

An illustration would be the room-temperature dual optical trap experiment of McCann, 
Dykman, and Golding ||, in which msi3x 10~ 16 kg and uj u = (7 ± 2) x 10 4 sec _1 . The 

corresponding force magnitude \jlmuj\kBT is approximately 10~ 13 Newtons. Any repetition 
of their experiment, with the addition of periodic driving, should take this dividing line into 
account. 

Another scaling-related issue has to do with the effects of choosing a period Tp for the 
forcing that is very small or large. In the Kramers limit, it is possible to take Tp oc e 13 , where 
j3 may be positive or negative. We make the natural choice (3 = 0, so that Tp is independent 
ofe. 

Preliminaries. — Our analysis of the a = case uses optimal trajectories. The e — > 
limit is governed by the action functional 

W [t ^ x(t)} = i J \x + V'(x) - fv(t)\ 2 alt. (4) 

First, suppose that / = 0. Then the most probable trajectory from x s to any specified 
point x' is the one that minimizes W [t \— > x(t)}. The minimum is taken over all trajectories 
from x s to x', and all transit times (infinite as well as finite). There is a single minimizer t h- > 
to each side of x s , which we term an optimal trajectory. The value W [t i— >■ x*(t)], which 
depends on the endpoint x' and may be denoted W(x'), is the rate at which fluctuations 
to x' are exponentially suppressed as e — » 0. In the steady state, the probability density 
p = p(x) of the particle will have the asymptotic form 

p(x) ~K(x)exp(-W(x)/e), e -> 0. (5) 

The prefactor K(x) must be computed by other means. 

Any such / = optimal trajectory must satisfy x = +V'(x), i.e., be a time-reversed 
relaxational trajectory. This is due to detailed balance, which holds in the absence of 
'hot spots'. The optimal trajectory from x s to x u is instanton-like: it emerges from x s at 
t = — oo and approaches x u as t — > +oo. Within the well, W(x) equals 2[V(x) — V(a; s )], so 
AE = W(x u ) equals 2 [V(x u ) — V(a; s )]. Also, K is independent of x. 

If / = 0, the model defined by the Langevin equation (fj) is invariant under time trans- 
lations. So the optimal trajectory from x s to x u is not unique. If x = x*(t) is a reference 
optimal trajectory, consider the family 
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t» x M{t)=x*(t + j-T F ), (6) 

Z7T 

where the phase shift 4> satisfies < <fi < 2tt, and t> is the period of the forcing function 
v = u(t). In the Kramers limit of any model with / nonzero but small, the most probable 
escape trajectory should resemble some trajectory of the form (|5|). That is, some m will 
be singled out as maximizing the chance of a particle being 'sloshed out'. A study of the 
/ — > limit should yield <f> m . 

This was the approach of ||. Suppose that / ^ 0. If AE is computed by applying 
(|) to the unperturbed (/ = 0) optimal trajectory x = x*(t), the first-order (i.e., O(f)) 
correction to AE will be fwi(<j>), where 

poo 

w 1 {<j>) = - xf\t)u{t)dt. (7) 



It is reasonable to average the Arrhenius factor exp (—AE/e) in the Kramers formula over 0, 
from to 2tt. If (•)</, denotes this averaging, then the escape rate will be modified by the 
driving, to leading order, by a factor (e~^' Wl ^ € )^. If a = 1, i.e., / = fie for some fi, then 
the Kramers formula (|3|) will be altered to 

f , M \/V"(x 8 ) \V"(x u )\ 

Clearly, m should be the phase that minimizes u>i (</>). 

Eq. (Bp is essentially the formula of Smelyanskiy et al. M. But our derivation makes it 
clear that their perturbative approach requires that / — > rapidly as e — > 0, i.e., that a be 
sufficiently large. Estimating the minimum of W[«] by applying it to unperturbed optimal 
trajectories yields a correction to AE which is valid only to 0{f 1 ). 

If / is independent of e, then by Laplace's method 

f e -f»iW*\ „ ■ 1 e i/2 e -fMM/* ( 9 ) 

as e — > 0. This would seemingly suggest that 

JV"(x s ) \V"(x u )\ l 
_ V — v wi 1/2 (—AE/e), (10) 

2vr J2*w>{(<i> m )f 



where A£? is shifted by fwi(<f) m ) to leading order, is the a = Kramers formula. But 
the prefactor in (|i~0|) is correct only in the small-/ limit. If / oc e a , the 0(f l ) correction 
to AE will be of magnitude e a . If a — 1, it will induce, as in ([Bp, a correction to the 
prefactor. But when a < 1/2, 0(f 2 ) corrections will also affect the prefactor. The most 
difficult case is a = 0, when computing the prefactor would require working to all orders 
in /. A nonperturbative treatment, like the analysis of Lehmann et al. [[|] or the following 
analysis, is needed. 

Analysis. — We first remove explicit time-dependence, when / ^ and Tp are fixed, by 
replacing (§) by 
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y = l. (11) 

Here < y < r F , and y is periodic: y = t f is identified with y = 0. The state space with 
coordinates X = (x, y) is effectively a cylinder. On this cylinder, the oscillating well bottom 
x = x^(t) is a stable limit cycle, and the oscillating well top x = x[f\t) is an unstable limit 
cycle. To stress /-dependence, we denote them Xy) and X^. 

To study escape through when e — > 0, we can employ results of Graham and 

Tel |TI|,[nj. The limit is governed by an instanton-like optimal trajectory X = X* (i) 
that spirals out of X^) and into X^l It is the most probable escape trajectory in the 

steady state. The exponent AE equals W[t i— > X* (t)], which in general must be computed 

( f) 

numerically. The trajectory X* would be computed nonperturbatively, by integrating 

Euler-Lagrange or Hamilton equations. 

( f) 

X* increasingly resembles a time-reversed relaxational trajectory, as it nears the oscil- 
lating well top. So at any specified y, the Ith winding of Xi , as it spirals into X^, has 
an inward offset n that shrinks geometrically, like ac~ l , as I — > oo. Here a = a(y) and c are 
/-dependent, and c = exp[/|V"(x<^(t))| dt\. 

The form (|5]) for the steady-state probability density generalizes 
to K(X) exp (— W(X.)/e). To compute W and K at any specified X', an optimal trajec- 
tory ending at X' is needed; in general, one different from xi^. An asymptotic analysis 
of the Smoluchowski equation for the probability density |5|,|13| shows that W satisfies the 
Hamilton-Jacobi equation 

(VW) ■ D ■ (VW)/2 + u • VW = 0, (12) 

and along any optimal trajectory, K satisfies 

K = -(V • u + D ij d i d j W/2)K. (13) 

Here u(x, y) = (—V'(x) + fv(y),l) is the drift on the cylinder, and (D^) = diag(l,0) is 
the diffusion tensor. It follows from (12") that the Hessian matrix (didjW) obeys a Riccati 



equation along any optimal trajectory |13|,|14[ . This gives a numerical scheme for computing 
K(K'). By convention, K is chosen to be Oil) on Xy). 

In principle, the steady-state escape rate A can be computed by the Kramers method [|TJ : 
evaluating the probability flux through X^. But this is intricate, due to a subtle problem 
discovered by Graham and Tel [P,^2]. Optimal trajectories that are perturbations of the 
escape trajectory t i— > X* (t) intersect one another wildly near X^. This is because 
t ^ xl /} (t) is a delicate object: a 'saddle connection' in the Hamiltonian dynamics sense. 
In consequence, any X' near X.\[> is reached by an infinite discrete set of optimal trajectories, 
indexed by I, the number of times a trajectory winds around the cylinder before reaching X'. 
The density asymptotics are |13| 



p(X) ~ $>«(X)exp(-H/«(X)/e), e - 0, (14) 
i 



since W and K are infinite-valued, not single-valued. 
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It is known [|TT|--|l3|j that at any fixed y, any is not quadratic but linear in the offset 



n 



from X^: 



W {l) (n) ps AE — |W jfOT | (aeT'n - (acT*) 2 ^) . (15) 

W }nn < is what, in the absence of multivaluedness, the Hessian matrix element d 2 W/dn 2 
would equal at n — 0. Along X^, it obeys the scalar Riccati equation 

dW, nn /dy = -W, nn 2 + 2V"(x ( u f) (y))W, nn . (16) 

W >fm = W >nn (y) is the 7>-periodic solution of this equation, which is easy to solve numeri- 
cally. At any y, W >nn equals 2V"(x u ) to leading order in /. Deviations from this value are 
due to anharmonicity of V at the well top. 

It is also known [jn| that the second term on the right-hand side of ( |13D tends rapidly 
to zero along X* , as it spirals into Xj/\ So with each turn, K is multiplied by exp[— /(V • 
u)dt], i.e., by exp[/ V"(x^\t)) dt]. This factor equals c -1 . So ~ Ac~ l for some A = 
A(y). Since n ~ ac~ l , it follows that along X* , K ~ k\n as n — > 0. Here &i = A/ a, 
like W^n, is a r^?-periodic function of y, which quantifies the linear falloff of K near X^. 
The linear falloff of K is a nonperturbative effect. 

As a function on [0, t f ), k\ turns out to be proportional to W^ nn [|1^]. It can be obtained 
numerically by integrating fll3|) along the trajectory ~X^/\ as it spirals into X^. It is the 
t — > oo limit of the quotient K/n. Deviations from constancy are due to anharmonicity 
of V. 

Substituting (|TJ) and ~ k x ac~ l into (0) yields 

oo 

e -AE/e j2 k iac - 1 exp {l^l [ac^n - (ac"') 2 /^ /e} (17) 

l=— oo 

as the e — > steady-state probability density p, at an inward offset n from the oscillating well 
top. Summing from — oo to oo is acceptable since the errors it introduces are exponentially 
small, and can be ignored. The dependence here on t, i.e., on y, is due to W >nn , k±, and a. 

Discussion. — The cycling phenomenon, and much else, follow from the infinite sum (p~7|) . 
To determine its behavior on the 0(e 1 ^ 2 ) diffusive length scale near the oscillating well top, 
set n = iVe 1 / 2 with iV fixed, and also multiply by e -1 / 2 . (As in the case of no periodic 
driving, a steady-state density p that is normalized to total probability 1 within the well 
must include an e -1 / 2 factor.) The resulting expression is invariant under e t— > c~ 2 e. So 

p(n = Ne 1/2 ,t) ~h i /\N,t)exp(-AE/e), e -> 0, (18) 

where the quantity h^(N, t), for any iV and any t in [0, Tp), is periodic in log e with period 
2 log c. 

In the steady state, the instantaneous escape rate X(t) through the oscillating well top, 
which equals (e/2)(d/dn)p\ n=0 , satisfies 

X(t) ~ (l/2)e 1 / 2 /ip(0, t) exp(-AE/e), e -> 0. (19) 

So at any t in [0,r F ), the instantaneous escape rate, divided by e 1//2 , ultimately oscillates 
periodically in loge with period 21ogc, i.e., with period 2[§\V" (x[f\t)) \ dt]. 
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Lehmann et al. noticed that on [0, 7>), the peak of the function A(«) may shift when 
e is decreased. Our results indicate that slow oscillations in the instantaneous escape rate 
are a widespread phenomenon. They have a simple physical cause. In the e — ^ limit, the 
most probable trajectory taken by an escaping particle is the helix 1 1— > X* (t), along which 
it moves in a ballistic, noise-driven way. However, once it gets within an 0(e 1 ^ 2 ) distance 
of the oscillating well top, it moves diffusively rather than ballistically. It is easily checked 
that the changeover to diffusive behavior takes place at a location that cycles slowly around 
[0, r F ), as e — ► 0. If e i — > c~ 2 e, the changeover returns to its original location. 

If the well top is perfectly harmonic, so that W >nn and k\ do not depend on t, and the 
bottom is too, it is straightforward to integrate X(t) over [0,t>). We find 



A ~ ^ .J,,/ - 6^exp(-A^/6). (20) 

y/2-KT F \V"(x u )\ 

It is useful to compare (|20| ) with the perturbative formula (JTQ|) . They can be reconciled if k\ 
diverges like /~ 1//2 as / — > 0. An f~ l l 2 divergence was seen in this special case by Lehmann 
et al. ||, and it occurs more widely ]T5| . It has major consequences. k\ is the normal 
derivative of the density prefactor K. But K is 0(1) on X^, and is well-behaved in the 
well interior as / — > 0. So there must be a layer near the well top, of width (9(/ 1 / 2 ), within 



which K slopes off to zero. The presence of this layer has been numerically confirmed |T5 

We can now compare the steady-state probability density flT8|), which is valid on the 
0(e 1 ^ 2 ) length scale near the oscillating well top, to the density when / = on the same 
length scale. The analog of h^\N,t), if / = 0, is 



erfc 



■y/\V"(x u )\N 



x e - 1/2 exp (\V"{x u )\N 



up to a constant. The first factor is the Kramers boundary layer function |J, and the second 
is from the Maxwell-Boltzmann distribution. 

It may seem odd that h^(N,t), which is defined by a complicated infinite sum, should 
degenerate into such a classical (and t-independent) form in the / — > limit. The details 
remain to be worked out, but the mechanism is clear: the / — > limit passes through an 
intermediate scaling regime, namely a — 1, where fll7l) does not apply. The dominant terms 
in the sum ( |PT|) are those for which ac~ l is comparable to e 1 / 2 . But in deriving fllTD , we 
used the linear falloff approximation: ~ k\ac~ l . As we saw, this is justified only if 

n = ac~ l J 1 / 2 . This will be the case for the dominant terms in the sum, provided that 
e 1 / 2 J 1 / 2 - So if a < 1, the formula flT7|) is valid in the Kramers limit. But if a > 1, it 
does not apply. 

In fact, the a = 1 case is a crossover regime, in which the 0(f 1 ^ 2 ) length scale is 
comparable to the 0{e 1 ^ 2 ) length scale. When / = f\e for fixed f%, the behavior of the 
0(e 1 ^ 2 ) boundary layer in the Kramers limit was determined by Smelyanskiy et al. ||. 
Presumably, their perturbatively derived expression interpolates between the boundary layer 
hy\», •) (as fi — ► oo) and the / = boundary layer of Kramers (as fi —> 0). 

In closing, we wish to emphasize the experimental importance of the scaling regimes with 
a > 0. Any system with periodic forcing / and noise strength e lies on an infinity of curves 
of the form / cx e a , indexed by a. It is the task of the experimenter to determine which of 
the corresponding Kramers limit behaviors, if any, applies. 

This research was supported in part by NSF grant PHY-9800979. 
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